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FOREWORD 


Contract  AF  04(6ll)-74U8  between  the  Air  Force  Flight  Test  Center,  Edwards  Air 
Force  Base,  California,  Air  Force  Systems  Command,  United  States  Air  Force  and  the 
Research  Laboratories  of  United  Aircraft  Corporation  pertains  to  Investigations  of 
the  characteristics  of  a  unique  gaseous-core  nuclear  rocket  concept.  These  investi¬ 
gations,  which  were  directed  toward  the  development  of  information  pertinent  to  the 
evaluation  of  the  feasibility  of  this  concept,  include  the  following;  (l)  experimental 
determination  of  the  characteristics  of  vortex  flow,  (2)  design  euid  fabrication  of 
equipment  to  provide  a  radial  temperature  gradient  in  a  vortex  flow,  (3)  theoretical 
calculations  of  the  characteristics  of  the  end-wall  boundary  layer  in  a  vortex  tube, 
and  (4)  theoretical  calculations  of  radiant  heat  transfer  in  the  proposed  engine. 

The  work  accomplished  under  the  contract  is  summarized  in  the  following  Research 
Laboratories'  reports  which  comprise  the  required  ? inal  Report 

I,  Summary  Report  (Report  R-2494-4) 

II.  Experimental  Investigation  of  Characteris  Ics  of  Confined  Jet-Driven 
Vortex  Flows  (Report  R-2494-2) 

III  Heat  Transfer  to  Confined  Vortex  Flows  by  Means  of  a  Radio -Frequency  Gas 
Discharge  (Report  R-2494-j) 

IV*  Theoretical  Solutions  for  the  Secondary  Flow  on  the  End  Wall  of  a  Vortex 
Tube  (Report  R-2494-1) 

V  Methods  for  Calculating  Radiant  Heat  Transfer  in  High-Temperature  Hydrogen 
Gas  (Report  M-1492-1,  Present  Report) 
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Report  M-1492-1 


Methods  for  Calculating  Realant 
Heat  Transfer  In  High -Temperature 
Hydrogen  Gas 

Contract  AF  OJ*(6ll)-744‘3 


SUMMARY 


'A  digital  compute:-  program  was  l-  velopei  to  detennine  heat  transfer  by  thermal 
radiation  in  any  gas  with  ^ne-dimensional  or  axisyrrmetric  temperature  and  pressure 
distributions.  Heat  transfer  war.  calcul  ;‘ei  for  a  one-dimensional  case  w'itli  two 
black-bj.iy  boundai-ies  and  a  m  inot  'nically  i  ucj'easing  tempf?i-ature  The  calculations 
u..  uitVi  Viv.i r"n£r»>»i  t  ^.runghout  the  sj ’lOi/  b'.'twoer.  the  Louijdaries  and,  '.vith  a  grey 

the  high-temperature  end. 


were  made  with  hydrogen  t.nroughout  the  sj 
gar.  added  to  the  hydrogen  a<- 

Tljir.  work  was 

ana  completed  as  part  cf  Cwntract 
Center,  |lU!'5b  Par.eT 

States  Air  Force. 


o'wn  research  program 

rest 

m 

ur.nnd,  United 


(.1)  In  gases  the  radiant  heat  flux  densities  temr.ine  1  by  a  diffusion 
analysis  may  bo  in  error  by  several  orders  of  magni  ule  for  carer  having  steep 
temperature  gradients  and  low  absorption  ooefficior.i  s . 

(2)  Accurate  vajucs  of  raaiant  heat  fi-ox  ue-nsity  can  bo  oDtained  for  all 
cases  by  employing  a  tranr,j'';rt  anaiy.'-.i.'.  ui.  i  uoing  ruffici- nt ly  small  increments 
in  the  numerical  integration  p.-'ucedure. 

('il  The  US',  of  a  liffusiun  analysis  f  or  caioul  ati. -nr  where  the  absorption 
coefficient  is  high  and  a  transport  analysis  where  r.he  abr.tvrpticn  coefficient  is 
low  is  t.he  most  practical  method  of  solution  for  a  general  digit. uL  computer  program 

The  b  rnid-free  absorption  of  the  H"  icn  is  1  l-.e  msrd  import -mi  -'ontri- 
butor  to  the  absorption  coefficient  of  hydrogen  ai  hX'  atm  bt'tw"'’r.  -  '  ’ )  K  and 
12,000  K. 


INTRODUCTION 


Quantitative  data  on  heat  transfer  by  thermal  radiation  are  required  In 
designing  high-performance  propulsion  systems,  electric  arcs,  and  hydrogen 
fusion  experiments  because  thermal  radiation  Is  an  Important  mechanism  of 
energy  transfer  at  high  temperatures.  Ihere  have  been  many  studies  pertaining 
to  radiant  heat  transfer  reported  in  the  literature  but  the  limitations  of 
each  studj'  preclude  their  use  in  systematic  calculations  for  many  practical 
cases.  For  example,  the  emisslvlty  of  hydrogen  'which  Is  of  particular  Interest 
in  propulsion  systems  due  to  its  low  molecular  weight)  has  been  calculated  In 
Ref.  1  for  temperatures  up  to  10,000  K,  and  the  continuous  emission  of  a  thin 
layer  of  hydrogen  for  a  range  of  temperatures  up  to  lb,600  K  has  been  calculated 
In  Ref.  2.  Althou^  the  latter  values  may  be  readily  converted  to  emisslvlty, 
this  quantity  Is  useful  only  In  caset-  Involving  layers  ■  f  gas  at  a  constant 
temperature. 

In  a  diffusion  analysis,  the  Rosselan'l  mean  opacity  and  the  temperature 
gradient  are  used  to  calculate  heat  flux  density.  The  Rosseland  mean  opacity 
of  hydrogen  has  been  calculated  In  Ref.  3  temperatures  between  6300  Bind 
50,400  K,  neglecting  the  Important  contributions  due  to  H"  bound-free  absorp¬ 
tion  and  H  free-free  absoxTiti^'n.  Unfortunately,  as  polntei  out  In  Chapter  V 
of  Ref.  4  and  proved  later  In  this  report,  a  diffusion  analysis  Is  not  appli¬ 
cable  to  gases  with  low  absorption  coefficients  and  steep  temperature  gradients. 

One-dlmenslonal  transport  analyses  of  thermal  radiation  in  a  gas  or  other 
semi-opaque  media  are  given  In  Refs.  5  and  6.  Neither  reference  contains  values 
for  hydrogen.  Reference  7  contains  a  one-dlmenslonal  transport  analysis  for 
pure  hydrogen.  However,  only  the  heat  flux  reaching  a  wall  Is  considered,  and 
the  extension  of  the  bound-free  H  continuum  to  lower  wave  numbers  due  to  the 
lowering  of  the  Ionization  potential  at  high  free-electron  densities  (Ref.  0) 

Is  neglected.  The  ’hange  In  composition  due  to  lowering  of  the  Ionization 
potential  Is  also  apparently  neglected,  and  a  less  accurate  cutoff  is  used  for 
H  states  than  In  Ref.  6.  /•  genersQ  transport  analysis  jf  thermal  radiation  In 

a  gas  is  given  In  Ref.  'll.  The  equations  developed  In  this  reference  have  not 
been  applied  to  hyirog'in,  and  their  solution  Is  much  more  complicated  than  the 
equations  for  a  one-dlmensional  transport  analysis. 

The  study  described  herein  was  Initiated  to  provide  a  generalized  digital 
computer  program  for  use  In  calculating  values  of  heat  flux  density  and  net 
energy  loss  for  any  gas  having,  one-dlmenslonal  or  axlsymmetrlc  temperature  and 
pressure  ilstrlbutlons.  In  developing  this  program  specialized  subroutines 
were  provided  to  expedite  calculations  for  hydrogen. 
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ANALYSIS 


Basic  Principles 

Thermal  radiation  is  composed  of  photons  moving  in  various  directions  with 
the  speed  of  light.  From  a  macroscopic  viewpoint,  it  is  convenient  to  define  a 
spectral  radiation  intensity,  Jgj  .  Consider  a  specific  direction  in  space  and 
an  infinitesimal  area,  ,  with  its  normal  at  an  angle  d  with  the  direction 
in  space  (Fig.  l).  Let  d^E  be  the  energy  of  the  thermal  radiation  traversing 
dA  in  a  solid  angle  about  the  specific  direction  in  space  C/Q.  ,  in  a  time 
d\  ,  and  in  a  wave  number  interval  dQ  .  Then  Jcj  is  defined  by  the  relation 
(Ref.  U) 


d^E  -  COS  9  d A  J £l  d\  dU) 


(1) 


so  that  Jq  ir  expressed  In  terms  -!;■  energy  per  unit  area  per  steradlan  per  unit 
time  per  unit  wave  number.  In  polar  c ordinates,  the  element  of  s-lld  angle, 

da  ,  is 


da  -  sin  (?  dd  d<^ 


(2) 


Consequently  the  net  energy  passing  through  dA  per  unit  area  per  unit  time  is 


0  0  0 

F  is  also  known  as  the  radiant  heat  flux  density. 

The  quantity  is  well  kr.rvn  for  a  black  body,  an  i  f .  r  thlr  case  Is 
represented  by  .  Fron  (Ref.  4) 


B<5 


2  h  ^ 
hcw 
-  HT 


(M 


Bw  is  Independent  of  the  angle  at  which  the  surface  of  the  black  L^dy  is  viewed. 
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For  a  gas  Jqj  is  more  difficult  to  determine  than  for  a  black  body.  Neglect* 
Ing  discrete  transitions,  there  are  three  basic  processes  which  affect  Jq  : 

(l)  two  particles  may  recombine  and  release  their  excess  energy  as  a  photon  In  a 
random  direction;  (2)  an  atom.  Ion,  or  molecule  may  absorb  a  photon;  and  (3)  two 
particles  In  a  field  of  photons  may  be  Induced  to  recombine  and  emit  a  photon  In 
the  direction  that  the  other  pliotons  sire  moving.  All  three  of  these  processes 
may  also  occur  when  a  free  electron  comes  close  to  an  atom.  Ion,  or  molecule.  In 
this  case  the  free  electron  either  loses  or  gains  kinetic  energy,  and  the  phenomenon 
Is  called  a  free-free  transition  or,  collectively,  Bremsstrahlung.  It  Is  convenient 
to  define  a  mean  free  path,  Xq  ,  for  radiation  of  wave  number  <jj  .  This  Is  the 
distance  that  would  be  required  In  a  uniform  gas  for  a  pencil  of  radiation  to  be 
reduced  to  jr  of  Its  original  value  If  only  absorption  is  considered  and  both  spon¬ 
taneous  and  Induced  emission  are  neglected.  The  reciprocal  of  Is  a  .  In¬ 

duced  emission  may  be  thought  of  as  negative  absorption,  so  an  effective  value  of 
a  ,  which  Is  designated  and  Is  less  than  a  ,  Is  (Ref.  A) 


hcG 


I-  e 


kT 


(5) 


It  Is  also  convenient  to  define  an  optical  depth  (Refs,  h  and  5)  as  follows: 


Tq  ■-  /qj  0's 


(6) 


Neglecting  spontaneous  emission  but  considering  absorption  and  Induced  emission, 
an  optical  depth  of  one  is  the  distance  required  for  a  pencil  of  radiation  to  be 
reduced  to  ^  of  its  Initial  value.  It  can  be  shown  fRef.  U)  that  J-  at  any 
point  In  a  gas  can  be  found  from  the  relation 


(7) 


which  Is  known  as  the  equation  of  tramsfer  and  Includes  the  effects  of  spontan¬ 
eous  emission,  absorption,  and  induced  emission.  Equation  (7)  applies  to  any 
specific  direction  in  space.  It  is  rigorously  valid  only  where  radiative 
equilibrium  exists,  but  its  use  elsewhere  is  merely  tantamount  to  the  customary 
assumption  that  the  three  Einstein  coefficients  are  independent  of  radiation 
intensity  or  flux  density  and  that  local  thermodynamic  equilibrium  exists. 


1+ 
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Diffusion  Analysis 

Equation  (7)  can  be  Integrated  directly.  The  value  of  Jq  for  T 


J 


w  ' 


Tc,—  -T 

/  Bjj  «>  “  flTTc  +C, 

-art 


where  C|  Is  a  constant  of  Integration.  It  Is  convenient  to  define  a 
optical  depth|  To'  ,  measured  baclcweurds  from  To.o  • 


Consequently 


d'i'  o'  +  c,  e 


-  T- 


S,  0 


But  the  location  of  the  axes  from  which  To/j  Is  measured  Is  arbitrary, 
cannot  depend  on  It.  Therefore 


-Tc  I 

0 


nay  be  expressed  as  a  Taylor  series. 


Bs  (Ts)  =  (0)-  T-  +.... 


Combining  Eqs.  (ll)  and  (12) 


(7B- 

Jo  =  (0) 


O' To’ 


3  le  then 

(8) 

new 

(9) 


(10) 

BO  Jo 


(11) 


(12) 


(13) 
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Note  that 


O' To’ 


o' Bp 
fl'S 


(IM 


"orablnlng  '13)  '!**) 


=  a^D  (0)  - 


dB 


(A) 


Ot 


ds 


■f  •  • 


(15) 


where  dS  io  aii  Inflnlteblnal  distance  In  the  same  direction  ab  Jjjj  .  Based 
on  Eq.  ;),  define  a  spectral  heat  flux  density  F-  In  the  x  direction 


,  2tr  .  w 


f'w.x  '•  Iq  Jui  Sin  B  cose  dd  di> 


(16) 


■’omblnlng  Eqs.  (1^))  and  (16) 


\  dB-  p 

F_  -Zrr  I  6-  (0)  sin  e  cose  de  ~2w  I  —  -7"“-  sin  0  cos  0  0^0+.... 

W|*  /(*)''  In  Ht 


Oj  c/x 


(17) 


4  w 


3a, 


dx 


"f  .  .  .  . 


(16) 


Pifferentlatlng  Eq.  (^^) 


dx 


2 


hCcii 

kT 


hci 


rfT 

dx 


(19) 
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Combining  Eqs.  (l8)  and  (19) 


/  -!!CC\^ 

3at  k  U-e  1 


(20) 


However,  the  first  term  In  Eq.  (20)  Is  not  always  sufficient  If  steep  tempera¬ 
ture  gradients  are  present  and  Is  very  small  In  certain  regions. 


One-Dlmensional  Transport  Analysis 


Rather  than 
exact  expression 
with  the  help  of 
a  normal  optical 


take  additional  terms  In  Bq.  (?0),  it  is  better  to  obtain  an 
for  the  radiant  heat  flux  density  for  the  one-dlraenslonal  case, 
Cq.  (7).  Insteai  of  measuring  optical  depth  In  all  directions, 
depth,  ,  Is  measured  adong  the  x  axis,  so 


T_ 


sec  B 


Combining  Eqs.  (7)  and  (21) 


sec  9 


(81) 


(82) 


The  solution  of  Eq.  (2i^)  Is 


'OJ 


'I 


To' 


-To"  seed  (  To'secd  „  -T o"  8  . 

J,-  =  g  /  sec  Be  “  +c,g  (83) 


u  '  '*2 


Considering  the  boundary  conditions  and  redefining  "  so  that  It  Is  measured 
backwards  from  the  point  of  Interest,  the  Intensity  In  an  Infinite  gas  Is 


sec  9 


(24) 
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Combining  Eqs.  (16)  and  (2U) 


r  ^  ■  -  /** 

F-  :  I  2Trsin0  COS0  /  B_  sec9~e 


-Tjj  "'sec® 


cfT.  dB' 

U) 


(25) 


T o  sec  9 ^ 


-  f  ^27rsin6'^  COS©'*'  f  sec  6*  9  °  a^T-  '  ^  dQ 

^0  •'o 


where  9  t  S  ^  t  T q  ,  and  are  shown  In  Klg.  2. 

Majtlng  the  substitution  y  ■  sec  9  and  char.glng  ‘■he  order  of  integration 


"  f 


>«o  fico  -T. " 


Ffij  :  2  TT  I  Bj^ 
'0 


"  y 


d^,dTa~-2v  I  I  dydTa'U26) 


'0  "'I 


Let 


yOO  _T-  '■ 

V 

E,  =  /  - ^ - dy 


(27) 


Combining  Eqs.  (26)  and  (27) 


--  2v  /*B5,  d  To'"-27r  1^"" 


(28) 


.alues  of  Eg  are  tabulated  In  Hef.  10  and  are  plotted  In  Elg.  j.  Physically 
Eq.  (26)  is  merely  the  difference  between  the  radiation  from  two  opposite 
hemisphere 6. 

In  order  to  make  Eq.  (28)  applicable  to  gases  extending  finite  dlstauices 
Instead  of  infinite  distances,  it  Is  necessary  to  consider  boundary  conditions. 
Assuming  the  boundaries  are  black  bodies  and  considering  absorption  smd  Induced 
emission,  the  radiation  at  the  reference  point  in  the  gas  due  to  Boundary  1 

Is  from  Eq.  (16) 
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R  OD  * 


/^2ir  sin  6~  cos  9~  e  dd 


(29) 


Substituting  y  «  sec  9 


>*’0  m 

-  y 

R.-.  ^  2ir  I  - -  flTy 


O) 


w 


I 


Let 


</y 


I 


(30) 


(31) 


Values  of  E3  are  tabulated  in  Ref.  10  and  plotted  in  Fig.  3»  Combining 
Eqs.  (30)  and  (3I) 


Rdi  =  E3 


(32) 


Considering  both  boundaries  and  combining  Eqs.  (26)  and  (32) 


’roj.'i 


T(d;'2 


'^u.  =  f  ’B^EjrfT."-  2-r  B.  Ej-z.  r  '  B^,  E ,  tf  Tj"*- 2.  B  E,  (33) 

^  r\ 


Evaluation  of  Eq.  (33)  requires  numerical  metnuac.  Equations  (17)  and  (18)  of 
Ref.  6  are  equivalent  to  Eq.  (33) 


Comparison  of  Diffusion  and  One-Dimensional 
Transport  Analyses 

Equation  (33)  always  gives  an  accurate  answer  if  the  numerical  analysis 
is  carried  out  with  sufficiently  small  Incrmnents.  Consequently  Eq.  (20)  can 
be  compeured  with  Eq.  (33)  to  obtain  an  idea  of  the  range  of  validity  of  the 
first  terra  of  Eq.  (20).  ^’o  do  this,  it  is  necessary  to  assume  the  relationship 

between  teraperature  and  normal  opticeJ.  depth.  For  a  linear  relationship,  the 
results  shown  in  Fig.  4  were  obtained.  The  temperature  Indicated  on  each  curve 
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was,  of  course,  the  temi)erature  at  the  point  for  which  the  flux  density  was 
calculated.  The  quantity  is  the  hemispherical  spectral  flux  density 

from  a  black  body.  It  is  evident  that  there  is  a  good  possibility  that  the 
first  term  of  Eq.  (20)  will  not  give  an  accurate  answer  if  F^j 

exceeds  0.3.  Unfortunately,  this  is  not  a  good  criterion  for  deciding  whether 
or  not  to  use  Eq.  (20)  or  (33)  because  the  criterion  is  sensitive  to  the  manner 
in  which  temperature  varies  with  normal  optical  depth.  In  certain  cases  the 
first  term  of  Eq.  (20)  may  give  answers  that  are  off  by  two  orders  of  magnitude. 


Axi symmetric  Transport  Analysis 

It  is  apparent  that  a  diffusion  analysis  is  not  adequate  for  all  cases 
likely  to  be  encountered,  while  a  one>dlmensionad  transport  analysis  is  based 
on  a  geometric  assumption  which  is  not  sufficiently  realistic  F-^r  some  propul¬ 
sion  systons.  A  rigorous  axisymmetrlc  transport  analysis  reqtilrec  solution  of 
the  relation 


F- 

<jj 


C  --sinScosac/e  d4> 


(34) 


which  is  the  r*'sult  of  combining  Eqs.  (11)  and  ^l6).  Solution  of  Eq.  (34)  is 
difficult  even  for  a  high-speed  digital  computer.  Instead,  it  appears  reason¬ 
able  to  Include  approximate  correction  factors  r/7  in  Eq.  (33)  to  allow  for 
the  decrease  in  flux  density  with  increasing  area  as  radius  increases.  The 
quantity  7  is  the  radius  to  which  Fq  applies  (also  the  reference  radius  for 
T  5j  "  )  and  is  sho-.Ti  in  fig.  Pecause  the  high-temperature  gas  near  the  center 
is  frequently  almost  opaque  and  has  only  moderate  temperature  gradients,  it  is 
frequently  a  good  approximation  to  replace  the  central  portion  of  the  gas  with 
a  black  body  inner  wall  at  the  same  temperature , as  shrwn  in  Fig.  With  the 
r  /  r  factors,  Eq.  (  3"^)  becomes. 


-2ir 


(35) 
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Equation  (3^)  glvas  a  reliable  result  for  heat  flux  reaching  the  outer  vail  If 
the  gas  Is  alaost  opaque  out  to  a  certain  radius  and  then  becoaes  transparent. 
This  Is  approximately  what  happens  In  typical  cases  with  hydrogen.  However, 
Hq*  (33)  should  be  used  with  great  caution,  as  It  does  not  give  zero  flux 
density  for  a  gas  at  unlfom  temperature  as  It  should  do.  Also,  It  gives  In- 
valid  answers  near  the  center. 


Over-all  Flux  and  Net  Energy  Loss 

To  find  a  sufficiently  accurate  over-all  radiant  heat  flux  density  (Includ¬ 
ing  radiation  of  all  wave  numbers).  It  is  necessary  to  integrate  the  spectral 
heat  flux  density  over  the  wave  ntsabers  where  thensal  radiation  could  possibly 
be  i^reciable  at  the  tenqMraturee  under  consideration. 


F  - 


d  o 


(36) 


nils  is  applicable  to  either  the  one -dimensional  or  axlsynnetrlc  case. 

Thr  .et  energy  loss,  €  ,  by  the  gas  due  to  spontaneous  emission,  induced 
emission,  and  absorption  of  thermal  radiation  of  all  wave  numbers  is 


div  F 
P 


(37) 


Consequently,  for  the  axi symmetric  case 


€ 


_l _ ^ 

r p  dr 


(rF) 


(38) 


The  Absorption  Coefficient 

To  detehnine  Cf  for  hydrogen,  it  is  necessary  to  specify  the  wave  number, 
temperature,  and  total  pressure,  and  to  calculate  the  partial  pressures  of 
ff,  H'*’,  H',  the  various  states  of  the  H  atom,  the  ionization  potentials  of  H“ 
and  H  ,  and  the  term  values  of  the  various  states  of  the  H  atom.  The  method 
for  determining  these  is  given  in  Ref.  8,  except  for  the  psurtlal  pressure  of  H" 
and  its  ionization  potential.  At  pressures  on  the  order  of  100  atm  and  temper¬ 
atures  at  or  above  6000  K,  It  is  permissable  to  assume 
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P  +  P. 


+  Ph* 


(39) 


and  solve  for  the  partial  pressures  of  Hg  ,  0  ,  h"*"  ,  and  H  as  specified  in 
Ref.  8.  P^-  is  lees  than  a  tenth  of  P^  and  less  than  4  x  10"^  p  ,  so  P^- 
can  be  found  to  a  good  approximation  by 


Ph  Ptf 
Ph- 


2  0Q^^  {2n  m)  (  k  T) 


HcIh- 

-TT" 


(40) 


where  1^-  Is  b.048  x  10^  cm*^  (Ref.  2).  The  variation  of  I  with  free- 
electron  density  was  not  considered  because  Ref.  1  states  that  the  lower 
effective  core  charge  of  the  H"  Ion  in  the  presence  of  free  electrons  pro¬ 
duces  a  smaller  change  In  Ionization  potential  than  for  an  atom.  At  100  atm 
and  temperatures  below  6000  K,  the  absorptioix  coefficient  of  hydrogen  does 
not  exceed  3  x  10*3  cra-^  below  109^000  wave  nmbers,  so  errors  due  to  the  use 
of  Eqs.  (39)  and  (40)  have  negligible  effect  on  heat  transfer  for  outer  walls 
up  t  ^  a  meter  In  diameter. 

There  are  nine  continuous  absorption  coefficients  for  the  species  and 
wave  numbers  Important  for  hydrogen  at  high  temperatures.  Five  of  these  are 
due  to  H  bound-free  transitions  where  H  has  a  principal  quantum  number  of 
from  1  to  5*  For  low  free-electron  densities,  the  relations  are  given  In 
Refs.  11  and  12.  However,  for  high  free-electron  densities,  the  Ionization 
potentled  of  H  Is  lowered  as  pointed  out  In  Ref.  8,  so  that  bound-free 
transitions  occur  at  smaller  wave  numbers.  The  equations  of  Refs.  11  and  12 
were  found  to  give  a  fair  approximation  for  these  smaller  wave  numbers.  Let 
a n  be  the  absorption  coefficient  (excluding  induced  emission)  for  bound-free 
absorption  of  H  with  principal  quantum  number,  n  ,  between  1  and  5*  Then 
If  (J  > 
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P  TT^  ff®  R 
H,n  y 


/3  k  T  h  ’  c  ^  w  ^ 


f  iTj 


1720  ,  —  -.1728 
\  Rv 


dRy  /w 


(W) 


If  d)  <  I  H  -  Wn 


On  = 


0 


(42) 
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Let  Qg  be  the  absorption  coefficient  (excluding  induced  emission)  for 
free-free  transitions  between  H'*'  and  a  free  electron.  Then  according  to 
Refs.  11  and  12 


4Ry 


3  ye  /?  T 


m 


5  5  j  5 

?  c^u,^  k2 


+  .1728 


kT  \3 


+  .3456 


i— 

hc<ii\R. 


(43) 


Let  a ,  be  the  absorption  coefficient  (excluding  Induced  emleslon)  for 
free-free  transitions  between  an  H  atom  and  a  free  electron.  Such  trsunsitions 
are  sometimes  called  H"  free-free  transitions.  The  pseudo  cross  sections  for 
such  events  have  been  evaluated  and  tabulated  In  Refs.  2  and  13  for  tempera¬ 
tures  between  2520  and  16,800  K  and  wave  numbers  between  548  and  2k, $90  cm"^. 
For  use  In  typical  cases,  these  .alues  were  extrajolated  to  k0,000  K  and 
27,420  cm"^.  "he  pseudo  cross  section  and  0^  are  related  by 


A 

H  H 


k  T 


(44) 


The  bound-free  absorption  coefficient  of  H“  Is  represented  by  Og  (ex¬ 
cluding  Induced  emission).  The  cross  section  for  this  transition  is  tabulated 
for  wave  numbers  between  546  auid  198,000  cra'^  In  Ref;^  I3  and  l4.  The  cross 
section  and  Og  are  related  by 


Ph-  ^  H- 
Pki 


(45) 


The  Hg  molecule  has  a  dissociation  c^jntlnuum  at  wave  numbers  greater  tlmui 
117,730  cra"^  'Ref.  15),  and  has  an  Ionization  contlnuur.  at  wave  nu-tbers  greater 
than  124,000  cm"4  (Ref,  16).  ""he  co>nblned  cross  section  In  megabarns  ( 10*1^  cm^) 
for  these  transitions  Is  given  In  Ref.  16  as  a  function  of  wave  length.  The 
cross  section  and  the  corresponding  absorption  coefficient,  Og  ,  are  related  by 


*9  = 


^2 

pkJ 


where  expressed  in  cm  . 


(46) 
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To  compute  a(  for  liydrogen^  the  induced  emission  factor  must  be  included. 
Neglecting  discrete  transitions 


0 


t  ' 


I 
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hco 

kT 


{^7) 


It  was  assumed  that  the  grey  gae  had  a  maximum  absorption  coefficient  of 
1  cm“^  for  all  wave  numbers,  rhe  particle  density  of  the  grey  gas  was  assumed 
to  be  small  compared  to  the  particle  density  of  hydrogen,  and  ionization  of 
the  grey  gas  was  neglected  In  computing  hydrogen  composition. 
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A  digital  computer  program  was  developed  to  use  either  (W)  or  (35)  to 
calculate  spectral  heat  flux  density.  Necessary  inputs  when  hydrogen  is  present 
include  ten^rature  and  pressure  at  eleven  e^uispaced  radii.  The  program 
utilizes  Eqs.  (39)  and  (UO)  together  with  the  procedure  described  in  Ref.  8  to 
calculate  camposition,  ionization  potential  cf  H  ,  and  term  values  for  H  .  The 
results  for  100  atm  and  12,C'0<)  :<  are  sunsnarlzed  In  Table  I.  The  absorption 
coefficients  for  hydrogen  are  then  calculated  by  Eqs.  fi*!'  through  (J»7)  end  are 
given  in  Fig.  6  for  100  atn  euid  12,0<'0  K.  Frovlblv^;-  is  hIso  made  to  add  arbitrary 
absorption  coefflclento  at  each  of  the  5  to  51  wave  numbers  used  at  the  eleven 
radii  to  take  into  account  grey  gas  or  other  gar.  ad  led  to  the  hydrogen. 

If  no  hydrogen  Is  present,  the  program  may  be  ured  for  any  gas  or  combination 
of  gases  by  computing  their  absorption  coefficients  independently  of  the  progrsm 
and  using  the  sura  of  these  absorption  coefflcientr,  ar  inputs  to  the  progiram.  In 
this  case,  T  must  be  given  at  eleven  radii,  but  P  is  net  needed. 

Next,  spectral  heat  flux  density  may  be  calculated  using  Sq.  (2C)  or  (35) 
as  applicable.  The  program  decides  whether  to  apply  Sq.  (20)  or  '35)  based  on 
the  magnitude  of  the  total  absorption  coefficient,  the  radius  increment,  and 
two  constants,  which  are  part  of  the  input.  Use  of  Eq.  I  35)  all  calcula¬ 
tions  would  unduly  complicate  the  program  and  increase  the  running  time  due  to 
the  large  number  of  integration  increments  required  if  Of  Is  high.  If  it  is 
desired  to  use  Eq.  (33)  instead  of  Eq.  (35)»  it  is  merely  necessary  to  add  a 
laurge  constant  to  all  eleven  radii.  Tlien,  In  Eq.  (  5) 
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80  Eq.  (35)  becomea  Eq.  (33)  for  all  practical  purposes,  and  a  one-dlnenslonal 
analysis  Is  obtained. 

When  all  the  spectral  heat  flux  densities  have  been  calculated,  Eq.  (36) 

Is  used  to  obtain  the  over>all  heat  flux  density.  Equation  (36)  can  then  be 
used  to  obtain  the  net  energy  loss  by  the  gas  at  eleven  radii. 

A  case  involving  hydrogen  and  a  grey  gas  was  calculated  baaed  on  the  teBqper- 
ature  distribution  given  in  Fig.  7,  a  pressure  of  1CX>  atm,  an  Inner  wall  radius 
of  I.32U  cm,  an  outer  wall  radius  of  15.2U  cm,  and  the  distribution  of  grey  gas 
indicated  by  the  absorption  coefficient  of  the  grey  gas  In  Fig.  8.  998*^76  cm 
was  arbitrarily  added  to  all  radii  to  obtain  a  one-dlmenslonsd  analysis. 

Equation  (33)  wau  usea  o  obtain  the  spectral  heat  flux  density  If  0^  was  less 
than  UU  cm'^,  whereas  Eq.  (20)  was  used  for  greater  than  cm~^. 

The  total  absorption  coefficient,  Qf  ,  Including  the  grey  gas  absorption 
ooefflclent,  is  shown  in  Fig.  8  for  a  wave  number  of  40,000  cm“^.  The  heat 
flux  density  is  shown  in  Fig.  9»  together  with  the  blaclcbody  hemispherical 
intensity  for  the  temperature  corresponding  to  the  radius.  The  blaek*body 
hemispherical  Intensity  is  the  heat  flux  density  that  a  black  body  would  radiate 
to  cold  surroundings.  The  heat  flux  density  exceeds  the  blaek>body  bemispberical 
intensity  at  radii  greater  than  about  8  cm.  The  beat  flux  density  reaching  the 
outer  wall  is  22  percent  of  the  inner-wall  black-body  hemispherical  intensity. 
Figure  10  shows  which  wave  numbers  are  important.  The  upper  curve  is  the  black- 
body  spectral  hemispherical  Intensity  for  the  temperature  of  the  inner  wall,  and 
the  lower  curve  is  the  spectral  heat  flux  density  reaching  the  outer  wall. 

Figure  11  gives  the  net  energy  loss  by  the  gas.  It  can  be  seen  that  the 
region  of  the  grey  gas  is  gaining  energy  from  the  hydrogen  at  small  radii,  whereas 
the  hydrogen  at  about  7  cm  Is  losing  energy  to  hydrogen  at  greater  radii  and  to 
the  outer  wall. 

A  second  sample  was  computed  without  adding  998.476  cm  to  any  radii.  The 
heat  flux  density  at  the  outer  wall  was  found  to  be  only  11  percent  of  the 
Inner  wall  black -body  hemispherical  intensity. 
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LIST  OP  SYMBOLS 

Absorption  coefficient,  l/cn 

2 

Area,  cm 

Black-body  intensity,  erg/em^  ster  sec 
Black-body  spectral  Intensity,  erg/ca  ster  sec 
Velocity  of  light^  cm/sec 
Constant  of  integration,  erg/em  ster  sec 
Constant  of  Integjratlon,  erg/cm  ster  sec 
Diameter,  cm 

Charge  of  electron,  esu;  also  base  of  natural  logarithms 
Energy,  ergs 

V 

A  definite  Integral,  dimensionless 
A  definite  integral,  dimensionleEB 
Radiant  heat  flux  density,  erg/on^  sec 
Spectral  racilcuit  heat  flux  density,  erg/co  sec 
Plancks  constsuit,  erg  sec 
Ionization  potential,  cm’^ 

Spectral  radiation  intensity,  erg/cm  ster  sec 
Boltzmann  constant,  erg/K 
Mass  of  electron,  gm 
Pressure,  atm 

Electronic  partition  function,  dimensionless 
Radius,  cm 
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7  R«ftr«ne«  radlu««  cm 

Ry  RydQMrg  coastMit  for  bgrdrogottf  em*^ 

Rq  Spectral  taemlepherleal  radiation  intensity,  erg/ca  sec 

$  Length,  cm 

t  Time,  sec 

T  Temperature,  K 

X  Distance  along  x  aids,  cm 

y  Sec  6  ,  dimensionless 

0  Conversion  factor,  *  lO*"^  atm  cm^/dyne 

<  Net  energy  lose  by  gas,  erg/gr  sec 

6  Angle  with  x  axis,  radians 

Spectral  mean  free  path  for  radiation,  cm 
p  Density,  gr/cm^ 

$■  Pseudo  cross  section,  cmV<*yne 

r  Cross  section,  ca^ 

Tq  Spectral  optical  depth  measured  from  fixed  axes,  dimensionless 

Tq'  Spectral  optical  depth  measured  backwards  frcmi  point  of  interest, 

dimensionless 

Spectral  normal  optical  depth,  dimensionless 
^  Azimuth  angle,  radians 

w  Term  value,  cm'^ 

S  Wave  nuBber,  cm*^ 

A  Solid  angle,  star 
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rterljts 

♦ 

Subscripts 

c 

OlFF 

9 

GRCY 

H 

m’ 

H2 

n 

0 

t 

TRANS 

V* 

X 

ui 


la  *  r  or  *  K  dlrsctlon 
In  -  r  or  -  K  direction 


Inner  wall 
Diffusion  analysis 
Free  electron 
Drey  gas 
Neutral  h 
Negative  h  Ion 
Positive  H  ion 
H2  molecule 

Principal  quantum  number  of  h 
Reference 

Including  Induced  emission 
Transport  analysis 
Outer  wall 
X  direction 
Spectral 
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TABLl  I 

Composition,  Ionisation  Potwitlal,  and  Ter*  Values 
for  Hydrogen  at  100  at*  and  12,000  K 


Pm 

”2 

•  3.3  X  10"^  at* 

Ph* 

■  I.5U  X  10"^  at* 

Ph* 

■  1.61  at* 

P# 

>  1.61  at* 

Pm,  ft  s  1 

>  9*6  X  10^  at* 

Ph,  ft » 2 

-  2.0  X  10*2  at* 

^H,  ft<S 

-  7.1  *  10’^  at* 

Ph,  ft*4 

-  5.5  X  10“*^  at* 

Ph,  ft  *  9 

•  0 

Ih 

•  1.0152  X  10^  cm’^ 

W; 

-  0 

-  8.226  X  10**  c*’^ 

w. 

•  9.772  X  10**  cm-i 

-  1.0412  X  10^  ca*^ 

"9 

-  1.1023  X  105  ca"^ 
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FIG.  2 


QUANTITIES  IN  A  ONE  >  DIMENSIONAL  TRANSPORT  ANALYSIS 


AND  Ex- DIMENSIONLESS 
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FIG.  3 


VALUES  FOR  Ee  AND  E3 
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ABSORPTION  COEFFICIENT,  CM” 
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ABSORPTION  COEFFICIENT  FOR  HYDROGEN 


TEMPERATURE 
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ASSUMED  TEMPERATURE  FOR  A  TYPICAL  CASE 

HYDROGEN  +  GREY  GAS 

P :  100  ATM 


0  :  30.48  CM  <  I  FT 
T(. :  19,444  K  :  35,000  R 
■  2222.2  K  :  4000  R 


0  2  4  6  8  10  12  14  16 

RADIUS,  r  -  CM 


BLACK  BODY  HEMISPHERICAL  INTENSITY,  ttB  -  ERG/CM^  SEC 
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HEAT  FLUX  DENSITY  FOR  A  TYPICAL  CASE 

HYDROGEN  -I-  GREY  GAS 


FIG.  9 


P  -•  100  ATM 
0  :  30.48  CM  :  I  FT 

Tc  =  19,444  K  s  35,000  R 
--  2222.2  K  :  4000  R 

ONE  DIMENSIONAL 


X  0  2  4  6  8  10  12  14  16 

RADIUE,  r  -  CM 
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spectral  heat  flux  density.  Fa,  OR  BLACK  BODY  HEMISPHERICAL  INTENSITY,  ttB*  -  ERG/CM  SEC 


FIG.  10 
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SPECTRAL  HEAT  FLUX  DENSITY  FOR  A  TYPICAL  CASE 

HYDROGEN  +  GREY  GAS 

P  «  100  ATM 


0  :  30.4S  CM  >  I  FT 
Tc>  19,444  K  >  35,000  R 
T^  «  2222.2  K  >  4000 R 
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NET  ENERGY  LOSS  FOR  A  TYPICAL  CASE 

HYDROGEN  4-  GREY  GAS 

Ps  100  ATM 
0  :  30.48  CM  :  I  FT 
Jq-  19.444  K  :  35,000  R 
T^  :  2222.2  K  :  4000  R 
ONE  DIMENSIONAL 


FIG. 


0  2  4  6  8  10  12  14  16 

RADIUS,  r-  CM 
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